Structure of hybrid protoneutron stars within the Nambu-Jona-Lasinio model 
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We investigate the structure of protoneutron stars (PNS) formed by hadronic and quark matter in /3- 
equilibrium described by appropriate equations of state (EOS). For the hadronic matter, we use a finite tem- 
perature EOS based on the Brueckner-Bethe-Goldstone many-body theory, with realistic two- and three-body 
forces. For the quarlc sector, we employ the Nambu-Jona-Lasinio model. We find that the allowed maximum 
masses are comprised in a narrow range around 1.8 solar masses, with a slight dependence on the temperature. 
Metastable hybrid protoneutron stars are not found. 

PACS numbers: 26.60.+C, 21.65.+f, 25.75.Nq, 12.39.-x 



I. INTRODUCTION 

Protoneutron stars are generally believed to be the evolu- 
tionary endpoints of the gravitational collapse of a massive 
star with mass larger than 8 solar masses |[ll|2]. Several differ- 
ent stages may happen during the evolution process Ini- 
tially, a PNS is very hot, with an entropy per baryon of the or- 
der of 1 to 2, and contains neutrinos produced by electron cap- 
ture in the beta-equilibrated matter. Due to their short mean 
free paths, they are prevented from leaving the star, and are 
temporarily trapped in it. The number of leptons per baryons 
that remain trapped is approximately 0.4. The subsequent evo- 
lution of the PNS, on a timescale of 10-20 s, is dominated by 
neutrino diffusion and cooling, and the newly formed neutron 
star (NS) stabilizes at practically zero temperature. 

In a previous article ^ we have studied static properties 
of PNS assuming that nucleons, hyperons, and leptons are 
present in stellar matter Our calculations are based on the 
EOS derived within the Brueckner-Bethe-Goldstone (BBG) 
theory of nuclear matter |6], extended to finite temperature. 
It turned out that for the heaviest PNS, close to the maxi- 
mum mass (about two solar masses), the central particle den- 
sity reaches values larger than 1 /fm^, thus making possible a 
hadron-quark phase transition. This could have strong conse- 
quences for the dynamical evolution of a PNS into a NS [7]. 

Therefore, we have extended the previous calculations, and 
studied hybrid protoneutron stars with a quark core described 
with the MIT bag model, finding maximum masses below 1 .6 
solar masses, independently on the temperature. Moreover, no 
metastable hybrid PNS were found [8]. 

In this paper, we discuss the structure of hybrid PNS us- 
ing the BBG EOS for describing the hadronic phase, and 
the Nambu-Jona-Lasinio model at finite temperature for the 
quark matter (QM) phase. We find that no phase transition 
to quark matter can take place if hyperons are included in the 
hadronic phase, with and without neutrino trapping. Thus, 
the hadronic phase contains only nucleons, and as a conse- 
quence, neutrino trapping will soften the equation of state, and 
the PNS masses will be smaller than the NS ones. 

We also find that the presence of QM limits the value of the 
maximum mass in a narrow interval around 1.8 Mq, with a 
slight dependence on the temperature. 



This paper is organized as follows. In section we 
review the baryonic EOS in the finite-temperature Brueckner- 
Hartree-Fock approach. Section |lll] concerns the QM EOS 
according to the Nambu-Jona-Lasinio model, whereas 
section |IV] contains results about the hadron-quark phase 
transition. In section [V] we present the results regarding PNS 
structure, obtained combining the baryonic and QM EOS for 
beta-stable nuclear matter Finally, in Section|Vl]we draw our 
conclusions. 



II. HADRONIC MATTER EQUATION OF STATE 
A. The BBG theory at finite temperature 

We start with the description of the hadronic phase. 
The EOS is based on the non-relativistic Brueckner-Bethe- 
Goldstone (BBG) many-body theory |6], which is a linked 
cluster expansion of the energy per nucleon of nuclear matter, 
well convergent [jgt] and accurate enough in the density range 
relevant for neutron stars. 

At finite temperature, the formalism which is closest to 
the BBG expansion, and actually reduces to it in the zero- 
temperature limit, is the one formulated by Bloch and De Do- 
minicis in |10]. This has been applied successfully to the 
study of the limiting temperature in nuclei |ll].In this ap- 
proach the essential ingredient is the two-body scattering ma- 
trix K, which, along with the single-particle potential U, sat- 
isfies the self-consistent equations 
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where ki generally denote momentum, spin, and isospin. W ~ 
Ej^^ +Eii^ represents the starting energy, = k^/2m + U (k) 
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the single-particle energy, and n{k) is a Fermi distribution. 
For the two-body interaction V, we choose the Argonne Vig 
nucleon-nucleon potential lfl2l] . We have also introduced 
three-body forces (TBF) among nucleons, adopting the phe- 
nomenological Urbana model [ 1 3.1 . This allows to reproduce 
correctly the nuclear matter saturation point po « 0.17 fm^^, 
E/A « —16 MeV, and gives values of incompressibility and 
symmetry energy at saturation compatible with those ex- 
tracted from phenomenology lfl4ll . For details, the reader is 
referred to Ref.fisl]. 

The calculation of stable configurations in PNS requires the 
knowledge of the EOS at different chemical compositions. In 
order to simplify the numerical procedure, we have introduced 
the so-called Frozen Correlations Approximation, i.e., the cor- 
relations at r 7^ are assumed to be essentially the same as at 
r = 0. This means that the single-particle potential Ui{k) for 
the component / can be approximated by the one calculated 
at r = 0. It has been shown in ref. that this assumption 
is valid with good accuracy if the temperature is not too high. 
Within this approximation, for given density and temperature, 
Eqs. ([TJ and ^ have to be solved self-consistently along with 
the equations for the densities, p,- — '^kni{k), and the free en- 
ergy density, which has the following simplified expression 



/ = E 



(3) 



where 



^/ = - L (ni{k)\nn,{k) + [\ - ni{k)]\n[\ - nt{k)]) (4) 
k 

is the entropy density for component ; treated as a free gas 
with spectrum Ei{k). In recent years, the BBG approach at 
zero temperature has been extended to the hyperonic sector in 
a fully self-consistent way lfl6l[l7l . by including the and 
A hyperons. For that, we have used the Nijmegen soft-core 
nucleon-hyperon (NH) potentials NSC89 |18], and neglected 
the hyperon-hyperon (HH) interactions, since so far no re- 
liable HH potentials are available. We have found that the 
presence of hyperons strongly softens the EOS, and produces 
a maximum NS mass that lies slightly below the canonical 
value of 1 .44 Mq ifTgll . However, since the quantitative effects 
of more reliable NH and HH potentials haven't been explored 
yet, in this paper we discuss finite temperature calculations 
by using non-interacting hyperons. This approximation could 
not be justified for the present work, because the hadron-quark 
phase transition occurs at densities much above the normal nu- 
clear matter saturation density, where hyperons are expected 
to play a role. A more complete study with the inclusion of 
interacting hyperons at finite temperature is in progress. 

For stars in which the strongly interacting particles are only 
baryons, the composition at given baryon density p is deter- 
mined by imposing electric charge neutrahty and equilibrium 
under the weak processes 
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where Bi and B2 are baryons and I is a lepton, either an elec- 
tron or a muon. When the neutrinos are trapped, these two 



requirements imply that the relations 
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are satisfied. In the expression above, x, = p,/p is the baryon 
fraction of the species /, bi the baryon number, and qj the elec- 
tric charge. Equivalent quantities are defined for the leptons 
I ~ e,pL. The initial PNS contains trapped neutrinos produced 
in electron capture process, so the electron and muon lepton 
numbers are conserved on dynamical time scales. We fix the 
electron lepton number = Xe + x^^ ~ 0.4, as indicated by 
gravitational collapse calculations, and i)i = — Xy^ = 0, 
since no muons are present when neutrinos become trapped. 

Hence, the composition of beta-stable and charge-neutral 
baryonic matter is determined by the baryon chemical poten- 
tials for each species /, which are related to the free energy 
density 



(8) 



Therefore, one needs to know the functional dependence of 
the free energy, Eq. ([3]), on the individual partial densities p, 
and on the temperature. In Ref. ISJ we have provided analyt- 
ical parametrizations of our numerical results for symmetric 
and neutron matter, from which one can readily obtain the 
nucleon chemical potentials in beta-stable matter The chem- 
ical potentials of the noninteracting hyperons and leptons are 
obtained by solving numerically the free Fermi gas model at 
finite temperature. 

Once the composition of the j3 -stable, charge neutral stel- 
lar matter is known, one can calculate the total free energy 
density / and then the pressure p through the usual thermody- 
namical relation 



•.dim 
dp 
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The resulting EOS is displayed in Fig. [T] where the pressure 
for beta-stable matter, without (solid lines) and with (dashed 
lines) neutrinos, is plotted as a function of the baryon den- 
sity at temperature T = 30 MeV. The upper curves represent 
the EOS for stellar matter containing nucleons only, whereas 
the lower curves display the case when free hyperons are in- 
cluded. We notice that if only nucleons and leptons are present 
in the neutrino-trapped matter, the EOS gets softer because 
nuclear matter is more symmetric. The presence of hyper- 
ons introduces additional softening. However, the degree of 
softening is smaller in the neutrino-trapped case than in the 
neutrino-free case, because the hyperons appear at larger den- 
sities in neutrino-trapped matter, and their concentration is 
smaller We have checked that this trend depends weakly on 
the temperature, which plays a role only in the low density 
range. 
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FIG. 1: The EOS in beta-stable baryon matter is displayed for tem- 
perature T = 30 MeV. The upper curves display results for matter 
containing nucleons and leptons, whereas the lower curves are for 
matter with hyperons. Solid (dashed) lines refer to the neutrino-free 
(neutrino-trapped) case. 



ni. QUARK MATTER EQUATION OF STATE 

A. The Nambu-Jona-Lasinio model 

We describe the quark phase within the Nambu-Jona- 
Lasinio model 120. 21,,]. The SU(3) version of the model in- 
cludes most of symmetries of the QCD Lagrangian, includ- 
ing the spontaneous breakdown of chiral symmetry, which 
is essential in treating the lightest hadrons. The NJL model 
assumes deconfined point-like quarks, and is not renormal- 
izable, requiring regularization through a cutoff in three mo- 
mentum space. In its commonly used form, the Lagrangian 
reads 

8 

- /i:{det/[^(l+75)^]+det/[^(l-75)^]}, (10) 

where q are the quark fields with three flavors and three col- 
ors, and Tfl are the U(3) flavor matrices. The model parame- 
ters are the current quark mass matrix m = diag(mj,,m£/,m.(), 
the coupling constants G and K, and the cutoff in the three- 
momentum space, A. Following ref.[20], we choose A = 
602.3 MeV, the two coupHng constants G = 1.835/A2 and 
K = 12.36/A^, as well as the current quark masses ot„ = = 
5.5 MeV, nis = 140.7 MeV. The values chosen for those pa- 
rameters have been adjusted to reproduce masses and decay 
constants of the pseudoscalar meson nonet. 

In the mean-field approximation at finite temperature and 
chemical potential, the thermodynamical potential of the 
quark phase is H = Hpc + ^int, where 



^(r,M,)=2M.r /^[ln(l-/,)+ln(l-/,)] 

^ i=u,d,s-' y'^^i 

(11) 

is the Fermi gas contribution arising from quarks. We consider 
three flavors and three colors, hence Nc ^3. The distribution 



functions of fermions and anti-fermions are given by fi = [l + 

exp(i(£',--^i))]"' and/i- = [1 +exp(i(£'; + ^i))]"'- Ei and 
/X; denote the single particle energy and chemical potential of 
the quark flavor /. 

The thermodynamical potential due to interactions among 
quarks is given by 



=11, d,: 
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2G < qiqi > 

4K < q^qu X qdqd >< q.sqs > 



(12) 



In the NJL model, the quark masses are dynamically gener- 
ated as solutions of the gap equation, obtained by imposing 
that the potential be stationary with respect to variations in 
the quark condensate < qtqi >, thus finding 



Mi 



- 4G < qiqi >^ +2K < qjqj > < qj^qk > (13) 



being {qi,qj,qii) = any permutation of (M,c/,i). The quark 
condensate < qiqt > and the quark number density n, are given 
respectively by 



< qtqi > = -2Nc 



= 2Nr 



d^p Mi 
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Eq.(fl4li has to be evaluated self-consistently with Eq.lflJTl. 
forming a set of six coupled equations for the constituent 
masses M,-. Once the self-consistent solutions are found, we 
can calculate the pressure and the energy density through the 
usual thermodynamical relations 



e = Q. + Ts + Y^^ini 



(15) 



where s — dQ./dT is the entropy density and = —d£l/dpLi 
are the number densities of flavor /. The total quark number 
density and the baryon number density are given by n = J^, 
and pB = n/3, respectively. The reader can find more details 
inref.f^l^. 

In a PNS with quark matter and trapped neutrinos we must 
impose beta equilibrium, charge neutrality, and baryon and 
lepton number conservation. More precisely, the individual 
quark chemical potentials are fixed by Eq. ^ with = 1/3, 
which implies 



(16) 



The charge neutrality condition and the total baryon number 
conservation read 



+ nfi = ^{2nu 
3PB = riu 



(17) 
(18) 



These equations determine the composition and the pressure 
of the quark matter phase. 
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Let us discuss first the main features of beta-equilibrated 
and charge neutral quark matter within the NJL model. In 
Fig. [2] we show particle fractions jc, = p,/p as a function of 
the baryon density for neutrino-free (upper and middle pan- 
els), and neutrino-trapped quark matter (lower panel). At 
zero temperature, we observe a substantial amount of u and 
d quarks, whereas the i-quark starts to appear at baryon den- 
sity four times larger than the normal saturation value. With 
increasing density, the s-quark population keeps smaller than 
the u and d populations. With increasing temperature, we 
observe a slight change of the particle population mainly at 
low density, because of the tails in the Fermi distributions. In 
this case, the onset for i-quarks takes place at density smaller 
than in the cold case. The presence of neutrinos influences 
quite strongly the composition, as displayed in the lower panel 
of Fig. |2] In this case the relative fraction of u quarks in- 
creases substantially from 33% to about 42%, compensating 
the charge of the electrons that are present at an average per- 
centage of 8% throughout the considered range of baryon den- 
sity, whereas the d and s quark fractions decrease. This im- 
plies that in the neutrino trapped case, being the strangeness 
content smaller than in the neutrino-free case, the equation of 
state for quark matter will be stiffer This is indeed shown in 
Fig. [3] where the pressure is displayed as a function of the 
mass-energy density for the neutrino-free (lower curves), and 
for neutrino-trapped case (upper curve). The behaviour of the 
pressure is consequence of the quark population in beta-stable 
matter. We clearly observe a kink in the pressure, correspond- 
ing to the onset of the i-quark. The kink becomes smoother 
in the neutrino-trapped case, because the .s-quark concentra- 
tion is smaller. Moreover, the EOS for /3 -stable, charge neu- 
tral quark matter derived within the NJL model shows a clear 
dependence on the trapped lepton fraction. This is different 
from the EOS derived with the MIT bag model for describing 
the quark phase, where no dependence on the trapped lepton 
fraction is actually found ll8l l22ll . 



IV. PHASE TRANSITION IN HOT BETA-STABLE 
MATTER 

In order to perform the hadron-quark phase transition in 
beta-stable matter at finite temperature, we adopt the sim- 
ple Maxwell construction. The more general Gibbs construc- 
tion ||23[ I24I] is still affected by many theoretical uncertainties 
ll25ll . and in any case the final mass-radius relation of massive 
(proto)neutron stars 1,26,1 is slightly affected. 

Assuming a first-order phase transition, we impose ther- 
mal, chemical, and mechanical equilibrium between the two 
phases. This implies that the phase coexistence is determined 
by a crossing point in the pressure vs. chemical potential 
plot, as shown in Fig. |4] There we display the pressure p 
as function of the baryon chemical potential /Xb for baryonic 
and quark matter phases at temperatures T = 0, 30 MeV. In 
the upper panel we show the case at T = without neutrinos. 
The solid line represents the calculations performed with the 
BBG method with only nucleons, the dot-dashed line displays 
the case when free hyperons are included, and the dashed line 
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FIG. 2: The fraction of u,d,s quarks is displayed as a function of 
baryon density for several cases, i.e. T = (upper panel), and T = 30 
MeV (middle panel) without neutrino trapping. The case for T = 30 
MeV with neutrino trapping is shown in the lower panel. 
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: Pressure as a function of the mass-energy density for beta- 
charge neutral quark matter. See text for details. 



is the NJL calculation. We observe that the phase transition, 
marked by a full dot, occurs if only nucleons are present in 
baryonic matter, as it was already found in ref. 12711 . In the 
BHF approach, this is due to the strong softening of the bary- 
onic equation of state when hyperons are included. The same 
holds true at finite temperature, without and with neutrino 
trapping, as clearly shown in Fig. ID In the middle panel. 
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FIG. 4: Pressure as a function of the baryon chemical potential for 
beta-stable, charge neutral nuclear matter and quark matter. See text 
for details. 



we display results at T = 30 MeV without neutrino trapping. 
We observe that the crossing point is only slightly affected 
by thermal effects, and is shifted towards smaller values of 
the chemical potential. This means that the phase transition 
starts at a baryon density smaller than in the cold case. Fi- 
nally, in the lower panel, we display the neutrino-trapped case 
at r = 30 MeV. Neutrino trapping makes the baryonic EOS 
softer and the quark matter EOS stiffer, hence the phase tran- 
sition is shifted at larger values of the chemical potential, and 
of the baryon density. 

In Fig.|5]we display the pressure as a function of the baryon 
density for the several cases discussed above. The plateaus are 
consequence of the Maxwell construction. Below the plateau, 
j3 -stable and charge neutral stellar matter is in the purely 
hadronic phase, whereas for density above the ones charac- 
terizing the plateau, the system is in the pure quark phase. We 
notice the neutrino trapping shifts the quark matter onset at 
larger values, and gives rise to a wider plateau. 

Similar calculations have been performed in ref.f28'l, where 
the non linear Walecka model [29,1 has been adopted for the 
hadronic phase, which includes nucleons and hyperons, and 
the NJL model with a different choice of the parameters for 
the description of the quark phase. The phase transition be- 
tween the two phases has been performed adopting the Gibbs 
construction. Also in this case, the onset of quark matter ap- 
pears at values of the baryon density which decrease with in- 
creasing the temperature. If neutrino trapping is taken into 
account, the onset of quark matter is shifted at larger values of 
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FIG. 5; Pressure as a function of the baryon density, normalized with 
respect to the nuclear matter saturation density po. 



the baryon density. In general, the larger the temperature, the 
smaller is the density of quark matter onset. We notice that, 
in spite of the different theoretical framework, we obtain quite 
similar results. 



V. PROTONEUTRON STAR STRUCTURE 

Based on these results for the beta-stable baryon and QM 
phases, we can now determine the properties of static (proto)- 
neutron stars. The stable configurations can be obtained from 
the well-known Tolman-Oppenheimer-VoUcoff (TOV) equa- 
tions IlU for the pressure p and the enclosed mass m, 



dp{r) 

dr 
dm{r) 

dr 



Gm{r)e{r) 



1 



£(,-) 



1 



n{r) 



1 ~2Gm{r)/r 



(19) 



once the EOS p{e) is specified, being £ the total energy den- 
sity (G is the the gravitational constant). Starting with a cen- 
tral energy density e(r = 0) = e^, the numerical integration of 
eqs.(fT9l) provides the mass-radius relation. 

Simulations of supernovae explosions fl isoll show that dur- 
ing the early stage, the entropy profile decreases from the sur- 
face to the core starting from values of 6-10 10], and the tem- 
perature drops rapidly to zero at the surface of the star due 
to the fast cooling of the outer part of the PNS, where the 
stellar matter is transparent to neutrinos. Therefore, we have 
modelled the PNS by assuming a hot isothermal and neutrino- 
opaque core separated from an outer cold crust lIsTl [32ll by 
an isentropic, beta-equilibrated, and neutrino-free intermedi- 
ate layer described by the Lattimer and Swesty equation of 
state OSll . For details, the reader is refeiTed to ref.i8l f34ll . 

In the following we consider two snapshots which repre- 
sent the thermodynamic al conditions in an evolving PNS : i) 
the initial state consisting of a PNS with a hot (Tcore ~ 30- 
40 MeV) neutrino-trapped core and a high-entropy transition 
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FIG. 6: The mass-central density (radius) relation is displayed in the left (right) panel. The mass is given in units of the solar mass Mq = 
1 .989 X lO-'-' g, and the central density Pc is normalized with respect to the saturation value po = 0.17 fm^-^ . 



layer (5env. ~ 6-8), joined to a cold outer crust; ii) the fi- 
nal state which represents the short-term cooling, where the 
neutrino-free core possesses a low temperature of about 10 
MeV and is direcly attached to a cold crust. 

The results are plotted in Fig.|6] where we display the grav- 
itational mass Mq (in units of the solar mass Mq) as a function 
of the central baryon density Pc (left panels) and the radius R 
(right panels). We display the complete set of results at core 
temperatures T = 0, 10 MeV without neutrino trapping, and 
T = 30, 40 MeV without and with neutrino trapping. Due to 
the Maxwell construction, the curves are not continuous and 
display plateaus in the mass-central density relation, which 
disappear in the case of the Gibbs construction. For values of 
central density smaller than the one characterizing the max- 
imum mass, the PNS are purely baryonic stars. Then an in- 
crease of central density is required in order to start the quark 
phase, as shown by the phase diagram of Fig|5] We find that 
the onset of the pure quark phase at the center of the PNS as 
the mass increases marks an instability of the star, i.e. the 
PNS collapses to a black hole at the transition point since the 
quark EoS is unable to sustain the increasing central pres- 
sure due to gravity. We had already found the same results 
in the case at T=0 IztIi . It has been argued that this instabil- 
ity might be related to the lack of confinement in the original 
NJL model [35.1 . Thus, heavy hybrid PNS in this model are 
practically baryonic stars up to large values of the central den- 
sities, eventually with a core in a mixed hadron-quark phase. 
The characteristics of the maximum mass configurations are 
reported separately in Table U We notice that the maximum 
PNS masses are comprised in a small range around 1.8 solar 
masses, with a slight dependence on the temperature. More- 
over, they turn out to be smaller than those of cold NS. This is 
due to the fact that the baryonic phase contains only nucleons, 
in which case neutrino trapping softens the EOS. This leads 
to smaller maximum masses, and excludes the possibility of 
metastable configurations. 



Composition 


T (MeV) 


Mc/Mq 


R(km) 


Pc/Pq 


£c{MeV/fm^) 







1.8 


10.3 


6.4 


129 A 


N,QM,l 


10 


1.82 


11. 


6.12 


nil 




30 


1.9 


15. 


6. 


1168 




40 


1.99 


23.8 


5.94 


1138 







1.76 


9.4 


7.82 


1751 


N,QM,l,Ve 


10 


1.78 


10.9 


7.76 


1730 




30 


1.86 


14.6 


7.7 


1689 




40 











TABLE I: Properties of the maximum mass configuration for differ- 
ent stellar compositions and temperatures. 

Below the maximum mass configuration, however, the stars 
develop an extended outer envelope of hot matter, the details 
of which depend on the treatment of the low-density bary- 
onic phase. The temperature dependence of the curves is 
quite pronounced for intermediate and low-mass stars, show- 
ing a strong increase of the minimum mass with temperature. 
Above core temperatures of about 40 MeV all stellar configu- 
rations become unstable. 



VI. CONCLUSIONS 

In this article we have studied hybrid PNS, combining the 
most recent microscopic baryonic EOS in the BBG approach 
at finite temperature with the Nambu-Jona-Lasinio model for 
describing the QM phase. Within this approach, we have 
found that the hadron-quark phase transition can take place 
if baryonic matter contains only nucleons. The inclusion 
of free hyperons softens considerably the EOS, and there- 
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fore inhibits any phase transition to quark matter. This is in 
agreement with the findings obtained at T = in Ref. llZTll . 
where the Njimegen soft-core nucleon-hyperon interaction 
was included. This requires new experimental data on hy- 
pernuclei, from which more modern parametrizations of the 
nucleon-hyperon and/or hyperon-hyperon interaction may be 
extracted. 

Within the NJL model for quark matter, we have also found 
the following : i) the transition density to the quark phase 
in neutrino-free baryonic matter occurs at about 6 times nor- 
mal nuclear matter saturation density; ii) with increasing tem- 
perature, the onset is shifted to smaller and smaller values of 
the baryon density, and iii) neutrino trapping substantially in- 
creases the transition density to QM. 

Therefore the hadron-quark phase transition studied with 
the original NJL model possesses general features, which are 
quite similar to the ones found with the MIT bag model ||^, 
where an increase of the temperature shifts the transition den- 
sity to lower values. Moreover, in this approach, hyperons are 
present in the hadron-quark mixed phase. 



We have also found that maximum mass for hybrid PNS 
lie in a narrow range around 1.8 solar masses, and that they 
are smaller than the corresponding ones obtained in the cold, 
neutrino-free case. Therefore, metastable configurations are 
not possible. This same result was found in ref.[8J, where the 
MIT bag model was employed for describing the quark phase. 
Both results are is in contrast to Ref. [f], where such metasta- 
bility was found for hybrid stars, but with a much stiffer bary- 
onic EOS including hyperons. 

Again, we confirm our prediction of limiting masses for 
PNS smaller than 2 solar masses, which are compatible with 
currently established observational data on NS. 



Acknowledgments 

We warmly thank J. Pons (University of Alicante, Spain) 
and M. Baldo (INFN Sezione di Catania, Italy) for fruitful 
discussions. 



[1] S. L. Shapiro and S. A. Teukolsky, Black Holes, White Dwarfs, 
and Neutron Stars (John Wiley and Sons, New York, 1983). 

[2] H. A. Bethe, Rev. Mod. Phys. 62, 801 (1990). 

[3] M. Prakash, I. Bombaci, M. Prakash, P. J. Ellis, J. M. Lattimer, 
and R. Knorren, Phys. Rep. 280, 1 (1997). 

[4] A. Burrows and J. M. Lattimer, Astrophys. J. 307, 178 (1986). 

[5] O. E. Nicotra, M. Baldo, G. F. Burgio, and H. J. Schulze, As- 
tron. Astrophys. 451, 213 (2006). 

[6] M. Baldo, Nuclear Methods and the Nuclear Equation of State, 
International Review of Nuclear Physics, Vol. 8 (World Scien- 
tific, Singapore, 1999). 

[7] M. Prakash, J. R. Cooke, and J. M. Lattimer, Phys. Rev. D52, 
661 (1995); A. W. Steiner, M. Prakash, and J. M. Lattimer, 
Phys. Lett. B486, 239 (2000); J. A. Pons, A. W. Steiner, M. 
Prakash, and J. M. Lattimer, Phys. Rev. Lett. 86, 5223 (2001). 

[8] O. E. Nicotra, M. Baldo, G. F. Burgio, and H. J. Schulze, Phys. 
Rev. D74, 123001 (2006). 

[9] B. D. Day, Phys. Rev. C24, 1203 (1981); H. Q. Song, M. Baldo, 
G. Giansiracusa, and U. Lombardo, Phys. Rev. Lett. 81, 1584 
(1998); M. Baldo, G. Giansiracusa, U. Lombardo, and H. Q. 
Song, Phys. Lett. B473, 1 (2000). 
[10] C. Bloch and C. De Dominicis, Nucl. Phys. 7, 459 (1958); 10, 
181,509 (1959). 

[11] M. Baldo and L. S. Ferreira, Phys. Rev. C59, 682 (1999); 
M. Baldo, L. S. Ferreira, and O. E. Nicotra, Phys. Rev. C69, 
034321 (2004). 

[12] R. B. Wiringa, V. G. J. Stoks, and R. Schiavilla, Phys. Rev. C51, 
38 (1995). 

[13] J. Carlson, V. R. Pandharipande, and R. B. Wiringa, Nucl. Phys. 

A401, 59 (1983); R. Schiavilla, V. R. Pandharipande, and R. B. 

Wiringa, Nucl. Phys. A449, 219 (1986). 
[14] W. D. Myers and W. J. Swiatecki, Nucl. Phys. A601, 141 

(1996); Phys. Rev. C57, 3020 (1998). 
[15] M. Baldo, I. Bombaci, and G. F. Burgio, Astron. Astrophys. 

328, 274 (1997); X. R. Zhou, G. R Burgio, U. Lombardo, H.-J. 

Schulze, and W. Zuo, Phys. Rev. 069, 018801 (2004). 
[16] H.-J. Schulze, A. Lejeune, J. Cugnon, M. Baldo, and U. Lom- 



bardo, Phys. Lett. B355, 21 (1995); H.-J. Schulze, M. Baldo, 
U. Lombardo, J. Cugnon, and A. Lejeune, Phys. Rev. C57, 704 
(1998). 

[17] M. Baldo, G. R Burgio, and H.-J. Schulze, Phys. Rev. C58, 

3688 (1998); Phys. Rev. C61, 055801 (2000). 
[18] P. M. M. Maessen, Th. A. Rijken, and J. J. de Swart, Phys. Rev. 

C40, 2226 (1989). 
[19] R. A. Hulse and J. H. Taylor, Astrophys. J. 195, L51 (1975); J. 

H. Taylor and J. M. Weisberg, Astrophys. J. 345, 434 (1989). 
[20] S. P Klevansky, Rev. Mod. Phys. 64, 649 (1992); P Rehberg, 

S. R Klevansky, J. Hufner, Phys. Rev. C53, 410 (1996). 
[21] M. Buballa, Phys. Rept. 407, 205 (2005). 
[22] D. P. Menezes, C. Providencia, and D. B. Melrose, J. Phys. G 

32, 1081 (2006). 

[23] N. K. Glendenning, Compact Stars, Nuclear Physics, Particle 
Physics, and General Relativity, 2nd ed. (Springer, New York, 
2000). 

[24] N. K. Glendenning, Phys. Rev. D46, 1274 (1992). 

[25] T. Endo, T. Maruyama, S. Chiba, and T. Tatsumi, Prog. Theor. 

Phys. 115, 337 (2006). 
[26] G. R Burgio, M. Baldo, P K. Sahu, A. B. Santra, and H.-J. 

Schulze, Phys. Lett. B526, 19 (2002); G. R Burgio, M. Baldo, 

P K. Sahu, and H.-J. Schulze, Phys. Rev. C66, 025802 (2002). 
[27] M. Baldo, M. Buballa, G. R. Burgio, R. Neumann, M. Oertel, 

and H.-J. Schulze, Phys. Lett. B562, 153 (2003). 
[28] D. P Menezes and C. Providencia, Phys. Rev. 068, 035804 

(2003); Phys. Rev. 069, 045801 (2004). 
[29] B. D. Serot and J. D. Walecka, Adv. Nucl. Phys. 16, 1 (1986). 
[30] J. A. Pons, S. Reddy, M. Prakash, J. M. Lattimer, and J. A. 

Miralles, Astrophys. J. 513, 780 (1999); L. Villain, J. A. Pons, 

P. Cerda-Duran, and E. Gourgoulhon, Astron. Astrophys. 418, 

283 (2004). 

[31] G. Baym, C. Pethick, and D. Sutherland, Astrophys. J. 170, 299 
(1971). 

[32] R. Reynman, R. Metropolis, and R. Teller, Phys. Rev. 075, 1561 
(1949). 

[33] J. M. Lattimer and R D. Swesty, Nucl. Phys. A535, 331 (1997). 



8 



[34] O. E. Nicotra, arXiv:nucl-th/0607055 pala', Phys. Rev. C75, 035804 (2007), and references therein. 

[35] M. Baldo, G.F. Burgio, P. Castorina, S. Plumari, and D. Zap- 



